Local 3D spin Hamiltonian as a thermally stable surface code 
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We study a 2D toric code embedded in a 3D Heisenberg ferromagnet in a brolcen-symmetry state 
at finite temperature. Stabilizer operators of the toric code are locally coupled to individual spins of 
the ferromagnet. The effects of the Goldstone modes of the ferromagnet in the ordered phase lead to 
an energy penalty for anyons that grows linearly with L, the linear size of the toric code. This 0{L) 
energy barrier for logical errors leads to a lifetime of the quantum memory that grows exponentially 
with L, assuming that the toric code is weakly coupled to a thermal bath with temperature below 
the phase transition of the ferromagnet. This provides a stable quantum memory with strictly local 
bounded-strength interactions in less than four dimensions. 
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Introduction. Topologically ordered phases of matter 
like Kitaev's toric code promise the possibility to store 
and process quantum information in a manner which is 
resilient against local imperfections [TH3]. However, a 
finite gap for the creation of topological defects (called 
anyons in the case of the toric code) is not enough to 
ensure stability against thermal fluctuations [Ujl]. If 
anyons can be created at a constant energy cost and prop- 
agate without any further energy penalty, they will at any 
non-zero temperature destroy the stored quantum infor- 
mation in a time which does not increase with the size of 
the memory. Indeed, it was shown that not only the toric 
code but a large class of 1-, 2-, and 3-dimensional Hamil- 
tonians suffer from the aforementioned thermal instabil- 
ity of quantum information [81110|. This is in contrast 
to the classical case, where magnetic devices allow the 
construction of self-correcting hard drives that are stable 
against both local perturbations and thermal excitations. 

Proposals for three-dimensional spin Hamiltonians 
with local few-spin interactions that do not fall victim 
to the aforementioned no-go results exist [TTHT^ . While 
the energy barrier in Refs. [12] and [T^ grows respectively 
like 0(log(X)) [H] and 0(L^/^), none of these models is 
expected to allow for a storage time increasing arbitrarily 
with system size. 

Several authors considered the possibility of suppress- 
ing the propagation or creation of toric code anyons by 
coupling them to external fields and found storage times 
that grow polynomially in L |15H19j or exponentially in 
L in a certain resonant parameter regime [T71 (THj . Sup- 
pressing the diffusion of anyons by means of an attractive 
interaction between them was proposed in Ref. and a 
possible implementation of such an interaction was found 
in Ref. [E] by coupling the toric code to a bosonic field. 
However, the coupling strengths in Ref. |15[ diverge poly- 
nomially in L. Note that quantum propagation of anyons 
due to magnetic perturbations in the Hamiltonian can be 
suppressed as well, namely through the presence of dis- 
order as shown in Refs. |20l|2Tj. Alternatively, Ref. [16] 
proposed long-range repulsion between anyons, leading 
to a vanishing anyon-density for large L. The implemen- 



tation and benefits of such long-range interactions have 
been further studied in Refs. [TTHTO] . 

In this work, we propose a three-dimensional spin 
model with purely local interactions of bounded strength 
that allows for thermally stable storage of quantum in- 
formation below some critical temperature. We consider 
a surface code embedded in a bulk ferromagnet and find 
an anyon chemical potential increasing linearly with L 
and thus a memory lifetime increasing exponentially with 
L. This scaling of the lifetime coincides with the four- 
dimensional toric code [H [H], which constitutes so far 
the only known example of a truly self-correcting quan- 
tum memory. The origin of this favorable behavior re- 
sides in the long-range character of the ferromagnetic 
transverse susceptibility in the ordered phase since it de- 
scribes the effective attractive interactions between sta- 
bilizer operators. This long-range behaviour is based on 
the Goldstone modes of the ferromagnet in the broken- 
symmetry phase. This result is general and does not 
depend on damping effects and temperature T as long as 
the ferromagnet stays in the ordered phase. While the 
long-distance behavior of the susceptibility can be un- 
derstood by rather general hydrodynamic arguments, we 
derive an explicit expression for the susceptibility in the 
one-magnon approximation for the sake of illustration. 
The effect of the ferromagnet can then be understood 
in terms of virtual emission and absorption of Goldstone 
modes (magnons) which mediate long-range interactions 
between the stabilizer operators. Finally, we show how 
to choose an appropriate coupling scheme to the ferro- 
magnet in order to hinder hopping of anyons (and not 
only their creation) by energy barriers that grow propor- 
tionally with L. 

Model. We introduce two sets of spins, namely Sj for 
the spins of the three-dimensional ferromagnet located at 
site j of a cubic lattice and Ij for the physical spins-1/2 
of the 2D surface code. The Hamiltonian we consider is 
purely local and takes the following form 

H = + Aj2WpS; , (1) 
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FIG. 1. A 2D surface code is embedded in the bulk of a 
3D ferromagnet which is in a broken-symmetry state. The 
stabihzers (plaquettes) Wp of the surface code locally couple 
to the S-spins of the ferromagnet (grey dots). Due to the 
broken symmetry, there exist long-range Goldstone modes or 
magnons (wiggly line) in the ferromagnet that mediate at- 
tractive gravitation-like interactions between the stabilizers. 



where A is the coupling constant between the spins of 
the surface code and the ferromagnet. Here, the pla- 
quette operator Wp = ^p_i-^p_2^p,3^p,4 is the poduct of 
spins around the square plaquette centered at p which 
are defined on a square lattice of linear size L with peri- 
odic boundary conditions (we set the lattice constant to 
unity). The three-dimensional vector p points towards 
the spin S lying in the center of a square plaquette of the 
surface code, see Fig. [Tj Note that this definition of Wp 
ensures that the blue and white plaquettes are equivalent 
to the usual toric code star and plaquette operators [1]. 

Hp = —JJ2{i j) Si -Sj with J > is the Hamiltonian of 
the three-dimensional Heisenberg ferromagnet with the 
sum restricted to nearest-neighbor sites The fer- 

romagnet is assumed to be in a broken-symmetry state 
pointing along the z-direction. In order to implement this 
formally, an infinitesimal magnetic field in z-direction is 
added to Hp. The physical spins Ij of the surface code 
are embedded in the bulk of the ferromagnet. The inter- 
action Hamiltonian Hint = ^X^p^p'^'p involves purely 
local five-body interaction terms. We define the anyon 



For A <^ J, we make use of a perturbative Schrieffer- 
Wolff transformation [23] similar to Ref. 24 to derive 
the effective plaquette-plaquette interaction (a detailed 
derivation is given in the Appendix) given by 
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where the coupling is Jp,p' = -A^Xxxip-p') and Xapir) 
is the static spin susceptibility of the S-spins only. Note 
that the energy in H^s can be minimized by either all 
stabilizers Wp having a -fl or a —1 eigenvalue. Adding 
the usual toric code Hamiltonian iJtoric = ~^X]p^p 
to Hamiltonian ([T]) would explicitly break this symme- 
try between anyons and anyon-holes (anyon-holes hp are 
defined through hp = 1 — rip). However, this term is 
not needed and its effect is vanishing in the limit of large 
L, as we discuss in the Appendix, so we neglect it for 
simplicity. The real space static susceptibility ^^^(r) is 
defined as the Fourier transform of 
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for a; = 0, where (...) denotes thermal equilibrium 
expectation values of the S-spins |25| at temperature 
T. The Fourier components are defined through = 
X- e~*'i'^'S'f , where Ng is the number of spins in 
the ferromagnet, and is a three-dimensional vector 
denoting the site of spin Si of the ferromagnet. 

It is not necessary to explicitly calculate the spin sus- 
ceptibility in the ferromagnetically ordered state to un- 
derstand its general behavior at large distances (or small 
q) pS]. Indeed, the spontaneous SO{3) symmetry break- 
ing of the state with finite magnetization pointing along, 
say, the z-axis, implies the presence of low-frequency 
Goldstone modes (called magnons in this context) and 
long-range correlations, i.e., the xx- (and yy-) suscep- 
tibility has to diverge for q and takes the follow- 
ing generic form in the hydrodynamic regime (low-energy 
and long wavelength regime) [26J 



Xxxiq,^^ = 0) 
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where i? > is the stiffness constant of the ferromag- 
net and = |(s^)p is the magnetization density with 
— ^ J2i Sf- The divergence at q —J' is trivially con- 
nected with the broken symmetry of the ground state: 
starting from a ferromagnetic state aligned along the 
z-direction, the slightest x-magnetic field is able to ro- 
tate and align all spins in x— direction and thus the re- 
sponse to an external magnetic field indeed diverges at 
q — 0. The same line of reasoning applies to the yy- 
susceptibility. Equation (|4| is thus very general and takes 
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into account damping effects due to magnon-magnon in- 
teractions |26j . Furthermore, an upper bound for the 
stiffness constant R can be derived [26], 
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Below we give an exphcit expression for R in the one- 
magnon approximation. 

The real space static susceptibility now follows by 
Fourier transforming Eq. Q which leads to 
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The plaquette-plaquette interaction is thus long-range 
gravitation-like and reads 
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Consequently, Eq. ^ describes a stabilizer Hamiltonian 
with long-range plaquette-plaquette interactions. Since 
the stiffness constant R is positive (see also below), the 
interaction between stabilizer operators Wp is attractive. 
By rewriting Hamiltonian ^ in terms of anyons we ob- 
tain up to some irrelevant constant (X]p p' -^p-p') 
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Note that for p = p' the two summands cancel, such 
that we can restrict the sums to p 7^ p'. A similar 
Hamiltonian has been derived in Ref. [I7j by coupling 
Kitaev's honeycomb model to cavity modes. While the 
non-locality of cavity modes allowed one to obtain con- 
stant long-range plaquette-plaquette interactions, here 
we start from a purely local Hamiltonian and obtain 
gravitation-like long-range interactions, i.e., Jp,p' c>c 

I — ^-TT- As we discuss below, this is sufficient to stabi- 
|p-p I ^ ' 

lizc the memory against thermal fluctuations. 

For the sake of illustration we calculate the stiffness 
perturbatively in the one-magnon approximation. For 
this we make us of the Holstein-Primakoff transforma- 
tion HZ] 
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in the limit rii ^ 25*, where — a^ai. Here, and 
al satisfy bosonic commutation relations and the associ- 
ated quasi-particles are called magnons. In Fourier space, 
we get Hp sa eqOqaq, up to some irrelevant constant, 
with magnon dispersion eq = 4JS'[3— (cos((72;) + cos(gj,)-|- 
cos{qz))], where we have introduced the Fourier com- 
ponents Gq = X^i 6~''''^'ai. Inserting Eq. M into 



Eq. ([3| leads to 
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Note that the temperature dependence drops out of the 
equation above since ([a±q,aj_q]) — 1. After a small- 
q expansion of eq (accurate if one is interested in the 
physics at large distances), we obtain 



xi°nq,^ = o) = 



2J|q|: 
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A comparison between Eqs. Q and ([TT| ) allows us to 
identify = 2JS^ since here (A/("^ = S^. Note 

that saturates the upper bound Eq. We thus 
obtain that Xxxi'r) = 8if7|r|' ^^'^ finally the following 
approximate plaquette-plaquette coupling 
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which is explicitly attractive since J > 0. We point out 



again that Eq. ( 12 ) is the one-magnon approximation of 
Eq. ([7]). Notably, the sole effect of both damping due 
to magnon-magnon interactions and temperature is to 



(8) renormalize the coefficients of the interaction (12 1, i.e 



(M(°))Vi?(o) ^ AP/R, while the form of the gravita- 
tional potential is not affected. Indeed, the diverging 
behavior of the transverse susceptibility at |q| — >■ 0, see 
Eq. Q, is a consequence of the broken S'0(3)-symmetry 
and is not modified by thermal fluctuations as long as the 
ferromagnet is in the ordered phase. Besides hydrody- 
namic and one-magnon approximations, other methods 
lead to results similar to Q for the spin susceptibility. 
In particular, the Tyablikov RPA-approximation 27 is 
especially useful for a spin- 1/2 ferromagnet and leads to 
the same q-dependence of the susceptibility independent 
of temperature in the ordered phase. We point out that 
the creation or annihilation of magnons will not gener- 
ate anyons in the toric code since the coupling terms in 
Hamiltonian ([ij commute with the toric code stabilizers. 
The presence of thermal magnons in the ferromagnet will 
only translate into small fluctuations of the chemical po- 
tential for the anyons. 

Note that the dimensionality of the ferromagnet used 
here is a central ingredient since Heisenberg ferromagnets 
in lower dimensions do not undergo a phase transition at 
T > ,28 . 

Thermally Stable Quantum Memory. In this section 
we show that the long-range attraction ([t]) between Wp 
operators leads to a stable quantum memory. Indeed, 
inserting Eq. ([t]) into Hamiltonian Q and using the fact 
that in the continuum limit X^p^o ~ Ilxl \p\ '^P ~ 
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cL, with a geometrical factor c € 0(1), we obtain 
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where the chemical potential of the anyons is 
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and thus increases linearly with the system size L. This 
implies a quantum information storage time that grows 
exponentially in L and /3, where /? = l/T is the inverse 
temperature of a bath weakly coupled to the memory, 
as originally demonstrated in Ref. [53] assuming that the 
interaction with the thermal bath can be described by the 
Davies equation. Below we present alternative arguments 
leading to the same conclusion. 



The second term in Eq. (13) describes a gravitation- 



like attraction between anyons, which is much weaker 
than the logarithmic potential found in Ref. [TS] (indeed, 
the energy needed for the anyons to escape the gravi- 
tational potential does not increase with L). Since this 
term helps to keep newly created anyon pairs attached 
to each other (for temperatures below the interaction 
strength A'^M'^ /ttR), it will have a further beneficial ef- 
fect on the memory lifetime. Note that this second term 
effectively reduces the anyon chemical potential. How- 
ever, this reduction is negligible since the anyon density 
is exponentially suppressed by the first term. 

Let us try to understand in more detail the deco- 
herence process of the memory in contact with a sim- 
ple model of a bath. We assume that the bath sup- 
ports single-spin processes in which an energy uj is trans- 
fered from the anyon system to the bath with rate 7(w) 
and that 7(0) 7^ [30]. Let S{N) denote the aver- 
age cost to create an anyon pair if there are already N 
pairs present. The gravitational interaction will lead to 
S{N > 1) < (5(0) = 2n{L) - A^1vP/[ttK). However, in 
the Appendix we show that this reduction will not lead to 
a finite self-consistent number of anyon pairs and that in 
fact we will have 5{N > 1) ~ (5(0) in the relevant regime. 

Since the presence of only two anyons diffusing across 
the memory leads to an uncorrectable logical error in 
times of order L^/7(0) [H], we need to show that the 
time for the creation of two nearby anyons that are not 
directly annihilated increases exponentially with system 
size. Whenever a new pair of anyons is created, their 
total hopping rate is given by 67(0) [3T] such that the 
probability that one of the two anyons ever moves be- 
fore the pair gets annihilated is 67(0)/[7((5(0)) -I- 67(0)]. 
Since 7((5(0)) = exp(/3(5(0))7(-5(G)) (which follows from 
the detailed balance condition) and the code consists of 
physical spins, we conclude that the total rate for cre- 
ation of anyon pairs that do not directly get annihilated 
is given by 

67(0) 



7(^(0)) + 67(0) 



< 6i^e" 



-I35(Q) 



7(0) . (15) 



The time needed to create such a pair is thus of or- 
der exp(/3(5(0))/L^7(0). In conclusion, we found a lower 
bound for the quantum memory storage time that in- 
creases exponentially with system size. 

This lifetime does not apply if the initial state of the 
system has anyons already present. Suppose that errors 
occur during preparation of the initial state, creating a 
finite density of anyons. If these errors are sufficiently 
sparse, it will be possible for error correction to recover 
the initial state. It is the job of the Hamiltonian to 
preserve this error correctability until the desired time 
of readout. The Hamiltonian given above will energet- 
ically favour the annihilation of anyons on neighboring 
plaquettes, undoing some of the errors. However, we can 
expect that a finite density of pairs will have been non- 
neighbouring, and so will remain. These only need to dif- 
fuse a constant distance to make correction ambiguous, 
which leads to a constant lifetime for the memory. To 
prevent this we can split the plaquettes into two types. 
'Strongly coupled' plaquettes are coupled to the ferro- 
magnet with a strength As . 'Weakly coupled' plaquettes 
have a strength A^ < Ag. These are chosen such that 
any sequence of single- or local two-spin errors that move 
an anyon from one weakly coupled plaquette to another 
must move it via a strongly coupled plaquette (exam- 
ple patterns are given in the Appendix). The chemical 
potential for the plaquettes will change from the form 



in Eq. (14 1, giving different values fisiL) and ybw{L) for 
the two types of plaquette. Performing the summation 
(as described in the Appendix) shows that the factor A^ 
becomes As A for HsiL) and A^^A for iiw{L) {A being a 
weighted average). The energy barrier required for anyon 
movement is therefore of order (1 — A^^/As)fJ.s{L), which 
increases linearly with system size. The resulting sup- 
pression of diffusion leads to a lifetime that increases ex- 
ponentially with system size, even when the initial state 
has a finite density of anyons. 

It may come as a surprise that associating some stabi- 
lizers with a lower energy penalty has a beneficial effect 
on the memory. However, note that the weakly coupled 
plaquettes allow energy to be dissipated from the anyons 
to the bath by hopping of an anyon from a strongly to 
a weakly coupled plaquette. On the other hand, if the 
chemical potential is independent of the anyon position. 



as in Eq. (14), this is only possible through annihilation 
of anyons. 

Conclusions. To the best of our knowledge, we have 
presented the first Hamiltonian in less than four dimen- 
sions that involves only local few-spin interactions of 
bounded strength and leads to a stable quantum mem- 
ory. We showed that the the presence of a ferromagnet 
stabilizes the memory by increasing the anyon's chemical 
potential linearly with the system size, and thus leads to 
a lifetime of the memory that grows exponentially with 
L. An anyon chemical potential that is proportional to 
L trivially implies a corresponding energy gap between 
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different groundstates of the toric code. This is not in 
contradiction to no-go results [8l-fl0] since our starting 
Hamiltonian ([T]) is not a stabilizer Hamiltonian and the 
effective stabilizer Hamiltonian ^ is highly non-local. 

The question of how to engineer the five-body interac- 
tions of Hamiltonian ([T]) remains open. One possibility 
may be to couple the spins of the ferromagnet to Kitaev's 
honeycomb model [32] (since it is equivalent to the toric 
code up to fourth order in perturbation theory), similar 
to Ref . [T7] . To conclude we would like to point out that 
our analysis is applicable to any stabilizer Hamiltonian 
and not restricted to the toric code. Indeed, the sole 
property which is needed to derive our results is that all 
stabilizer operators Wp commute among themselves and 
get coupled by massless Goldstone modes. 
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APPENDIX 

Interactions mediated by a translationally invariant 
system. For the sake of completeness, we show here a 
detailed derivation of Eq. (2) of the main text with the 
use of a perturbative Schrieffer- Wolff transformation sim- 
ilar to Ref. [T]. Let us thus consider the Hamiltonian 

H ^Ho + V = Ho+J2 "^q^-q > (16) 
q 

where Hq is a general S-spin Hamiltonian and Ai arbi- 
trary operators which commute with Hq and with each 
other. The Fourier components are defined through 
^q = 7^E.e-''i-^'S, and A^ = ^J^^e-^'^^A, 
where denotes the number of spins and their 
site. The effective Hamiltonian up to second order in A 
is then simply given by 

H^^ =Ho + U = Ho--lun dt e""* [V{t) , V] , 

(17) 

where V{t) = e'^"*t^e '^"t jg given in the Heisenberg 
representation (h = 1). If we now insert the explicit 
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expression for V we obtain 

POO 

q,q' 

= "2 ^E / e-"*([^^(i), 

q.q' ° 

+ .(18) 

=0 

We assume that the S-spins are in thermal equilib- 
rium, described by the canonical density matrix p — 
^-0Hf /Tr e"*^^^ , where Hp is the S-spin Hamiltonian 
without the coupUng to the plaquettes and corresponds 



to the main part of the Hamiltonian in Eq. (16), i.e.. 
Hp — Hq. In doing so, we neglect the backaction of the 
surface code on the ferromagnet. This is justified for two 
reasons. First, to include the coupling term on A in p 
would lead to higher order corrections in the coupling 
constant A (we retain only the leading order terms up to 
order A'^ in the effective Hamiltonian Hcs). Second, the 
surface code, being a two-dimensional lattice, acts only 
as a boundary perturbation on the three-dimensional fer- 
romagnet, and thus this perturbation becomes negligibly 
small in the thermodynamic limit. Further, formally we 
need to add an infinitesimal fictitious Zeeman field to Hp 
that breaks the 50(3) symmetry explicitly and induces a 
spontaneous ferromagnetic order below a critical temper- 
ature. Such a field also regularizes the zero-mode of the 
ferromagnet and thus makes the Schrieffer- Wolff series 
well-defined formally. After taking the thermodynamic 
limit, we then let this field go to zero in expectation val- 
ues of physical quantities. Here, we rely on a formal 
perturbation expansion in powers of A/ J. Convergence 
of this formal expansion is an interesting question by it- 
self and can be approached along the lines lines discussed 
in Ref. [5]. However, such rigorous treatment is beyond 
the present scope and will be addressed elsewhere. 

The equilibrium expectation values are then denoted 
by (...). If Hq is translationally invariant, such 
that (S^^S^^) = {SgS^^_^J, we have (S^S^,) = 
(S'qS'"q)(5q+q',o, and thus 



u = 



hm ^ / dt e-"* {[S^it), S%]) AqA_q 



(19) 



where Xaai'i) is the static spin susceptibility. 

Mean-field analysis of anyon pair creation. Assume 
that there are already N anyon pairs present. We want 
to determine the average (averaged over all possible po- 
sitions of the existing anyons) energy cost S{N) to create 
a new pair. From the point of view of one of the two 
newly created anyons, we assume that the existing 2N 



anyons are uniformly distributed over all — 2 remain- 
ing positions. The averaged interaction between one of 
the newly created anyons and each existing one is thus 



X2 



^— (^■2'£j^^o + AH'iy{7rR)^ 



L2-2 



{2piL)-AHiy{nR)) , (20) 



where we have subtract the energy — Al'^ / {tt R) due to 
attraction with the other anyon of the same pair. In- 
deed, we are only intersted in the attraction energy due 
to anyons which are already present before the creation 
of the pair. The total energy S{N) to create the new pair 
is thus given by 



6iN) 



SiO)- 

m 



4N 

AN 



(2^(L) - A^M'^/i^K)) 



1 



i2 



(21) 



where (5(0) = 2p{L) - A'^AP/{ttR). 

The mean-field energy of N anyon pairs is thus 



6{0)N- 



2N 



(22) 



The symmetry N /2 — N is, reminiscent of the fact 
that the energy in Eq. (2) of the main text can be min- 
imized by either all stabilizers having a -|-1 eigenvalue 
(no anyons present) or a —1 eigenvalue (memory full of 
anyons). The energetic gap between the sector in which 
there are almost no anyons and the sector in which the 
memory is full of anyons is of order A(0)L2 = 0{L^), 
so transitions between these two sectors happen on time- 
scales much longer than the time before the stored quan- 
tum information is lost. Consequently, each sector may 
serve as a thermally stable quantum memory, but at each 
moment in time we can only use one of the two. Without 
loss of generality, we consider the case where the sector 
with (almost) no anyons present is used for quantum in- 
formation storage. 



From Eq. ([21]) we have that 5{N) = (5(0)(l-2n), where 
n denotes the density of anyons. As there can only be 
zero or one anyon at each position, we obtain the self- 
consistent equation for the mean-field anyon density in 
equilibrium, 



[exp (/3(5(0)(l-2n„,f)) + l]- 



(23) 



If the left-hand side of this equation is smaller/larger 
than the right-hand side, the anyon density will tend to 
increase/decrease. If rimt solves this equation, so does 
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1 — n^f. One self-consistent density is Umf = 5- The sta- 
bility of this density depends on the temperature of the 
bath. For fi6{0) < 2 we have a unique self-consistent den- 
sity rimf = i and this density is also stable. For (3S{0) > 2 
the density ^ becomes unstable and two new stable self- 
consistent densities n* and 1 — n* emerge (let n* denote 
the smaller of the two) . The system of anyons therefore 
shows a phase transition and spontaneous breaking of 
the anyon anyon-hole symmetry at a critical temperature 
(5(0)/2, which is of order ^L. However, this tempera- 
ture becomes only relevant if it is lower than the criti- 
cal temperature of the ferromagnet which is of order J, 
which will not be the case in the limit of large L. For the 
purpose of quantum information storage, we are clearly 
interested in temperatures below both of these critical 
temperatures. 

Adding the toric code Hamiltonian i?toric ~ 
— A^p Wp to Eq. (1) in the main text explicitly breaks 
the symmetry between anyons and anyon holes and will 



lead to an additional summand 47V A in Eq. (22). How- 
ever, the modification of the self-consistent densities n* , 
1 — n* , and ^ through this new term becomes vanishing 
for large L, as A does unlike (5(0) not grow with L. 

Let us consider the self-consistent solution n* . We 
want to show that n* is exponentially suppressed with L 
and consequently that the number of anyons itself goes to 
zero in the thermodynamic limit. After straightforward 
algebra, one can show that n — 2e^^^^^^ < 1/2 with 



/3(5(0)e- 



PSiO) ^ log(2) 



(note that this condition is readily 



satisfied since S{0) grows linearly with L) satisfies 
[exp(/3A(0)(l - 2n)) + 1]-^ < n , 



(24) 



and therefore n > n* . Since n is by definition exponen- 
tially suppressed with L and n* < n we finally conclude 
that the self-consistent solution n* of Eq. ( 23 1 goes ex- 



ponentially to zero with L. A direct consequence of this 
is that the equilibrium number of anyons n* also van- 
ishes exponentially with L and will generally be much 
smaller than the minimal positive value 2. Hence the 
anyon number will fluctuate between and small even 
integers, such that S{N) « 5(0) from Eq. (pTi. 



An example pattern for strongly and weakly coupled 
plaquettes. In the toric code model there are two types 
of anyon, e and m, which reside on two kinds of pla- 
quette, s and p, respectively. Note that, when the code 
is defined with spins on the edges of the lattices, these 
correspond to the stars and plaquettes, respectively. 

Consider a spin in the square tiling of Fig. [2] (a), shared 
by two s-plaquettes to the top-left and bottom-right and 
two p-plaquettes to the top-right and bottom- left. The 
application of a Pauli to such a spin will affect the e 
anyon occupations of the two s-plaquettes. If both were 
initially empty, an anyon pair will be created. If both 
initially held an anyon, this pair will be annihilated. If 
only one held an anyon, it will be moved to the other 




FIG. 2. Two tilings of plaquettes are shown on which the 
code may be defined. Spins are located on vertices, (a) The 
square tiling, as used in the main text, s-plaquettes are shown 
in dark blue, blue, or light blue, p-plaquettes are shown in 
white, (b) An alternative tiling, with alternating triangular 
and octagonal plaquettes. s-plaquettes are shown in dark blue 
and blue, p-plaquettes are shown in white and grey. 



plaquette. The application of a Pauli 1^ has the same 
effect for the m anyons of the p-plaquettes. For spins 
where the positions of s- and p-plaquettes are exchanged, 
the roles of I" and P are also exchanged. No operation 
exists that can move an anyon from an s-plaquette to a 
p-plaquette, or vice-versa. 

Creation, movement and annilation of anyons are 
therefore achieved by Pauli operations. Using single spin 
operations, creation of a pair will always lead to the 
anyons occupying neighboring plaquettes (where neighor- 
ing means that they share exactly one spin). Similarly, 
single spin operations can only move anyons from one pla- 
quette to a neigboring one, or annihilate anyons on neigh- 
boring plaquettes. Since we assume that the system- 
bath coupling supports only single spin errors, it is ex- 
actly these processes that we consider during thermaliza- 
tion. However, it should be remembered that two-spin 
perturbations may also be present in the Hamiltonian. 
Local two-spin errors should also be considered, which 
can create, annihilate and transport anyons on next-to- 
neighboring plaquettes. 

With this in mind, we wish to split both s- and p- 
plaquettes into two groups, one of which will be strongly 
coupled to the ferromagnet with a coupling and the 
other of which will be weakly coupled with a strength 
A^, < Ag. This will give the plaquettes of the former a 
higher chemical potential than those of the latter, with 
an energy difference that increases linearly with system 
size. 

The pattern of strongly and weakly coupled plaque- 
ttes should be chosen such that anyons become trapped 
within the latter, which will occur if two conditions are 
satisfied. Firstly, any anyons initially on strongly cou- 
pled plaquettes should quickly move into a nearby weakly 
coupled plaquette. Secondly, it should not be possible 
for anyons to be moved from one weakly coupled pla- 
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quette (or a small cluster of weakly coupled plaquettes) 
to another by a sequence of either single- or two-spin 
operations without passing through a strongly coupled 
plaquette. 

The first condition can be met if anyons on strongly 
coupled plaquettes cannot be moved large distances by 
a sequence of either single or two spin operations with- 
out either moving through a weakly coupled plaquette, 
or through a strongly coupled plaquette that neigbors a 
weakly coupled one. The latter is relevant because the 
difference in chemical potential between the two types 
of plaquette will exponentially suppress the probability 
that the anyon will do anything but decay into the weakly 
coupled plaquette in this case. 

Both conditions are satisfied by the pattern shown in 
Fig. [2] (a). Here, weakly coupled s-plaquettes are shown 
in dark blue. Strongly coupled s-plaquettes that neigh- 
bor weakly coupled s-plaquettes are shown in blue, and 
those that do not are shown in light blue. Regions of 
strongly coupled plaquettes that do not neighbor weakly 
coupled plaquettes are separated from each other by a 
width of three spins. Sequences of one- and two-spin op- 
erations therefore cannot move anyons in one such region 
to another without going via strongly coupled plaque- 
ttes that do neighbor weakly coupled plaquettes, which 
will almost certainly result in the anyon decaying into 
the neighboring weakly coupled plaquettes. Similarly, re- 
gions of weakly coupled plaquettes are separated by the 
same width, preventing movement between them without 
going via strongly coupled plaquettes. 

The initial movement of anyons on strongly coupled 
plaquettes to nearby weakly coupled plaquettes may 
cause ambiguity for error correction if the error rate dur- 
ing initialization is too high. Even so, for sufhciently 
low error rates this movement will have no effect on cor- 
rectability. Once the movement is complete, the expo- 
nential suppression of diffusion will then ensure that the 
correct ability of the errors is preserved for a time expo- 
nential with the system size, since such an exponentially 
long timescale will be required for the anyons to climb 
out of the weakly coupled plaquettes. 

We will now demonstrate that the difference in chemi- 
cal potentials between strongly and weakly coupled pla- 
quettes leads to the energy barrier required to suppress 
diffusion. To determine the chemical potential of an ar- 
bitrary plaquette p (which is either s- or p-type), the 
following sum over all plaquettes must be performed 



takes the form 



Mp(i) 



2^ 



5Z 

p' 



P'l 



(25) 



where the prime in ^ means that p' 7^ p. Here A^i 
denotes the coupling of plaquette p' which will be As 
or depending on whether this plaquette is weakly or 
strongly coupled, respectively. By numerically perform- 
ing the summation we find that, in the L — ^ 00 limit, it 



p' 



A. 



3A, 



p' " 



P'l 



c-L 



(26) 



where c is an 0(1) constant, as mentioned in the main 
text. Its value does not depend on whether the sum is 
centered on a strongly or weakly coupled plaquette, or 
at least does not do so to a non-negligible degree. To 
three decimal places its value is found numerically to be 
c = 3.524. The linear combination of Ag and A^ is a 
weighted average A — {iAg -f A^)/^, which arises from 
the fact that there are three times as many strongly cou- 
pled plaquettes as weakly coupled plaquettes. The chem- 
ical potentials for weakly and strongly coupled plaquettes 
are then 



cA.AM^ 



f^w{L)^ -L. (27) 



Clearly, fJ-s{L) — fJ-wiL) = 0{L), giving the required en- 
ergy barrier. 

Alternative tiling with four-body coupling. A pattern 
of strongly and weakly coupled plaquettes, stable against 
single-spin errors, is shown for an alternative tiling in 
Fig. [2] (b). Strongly (weakly) coupled s-plaquettes are 
shown in blue (dark blue) and strongly (weakly) coupled 
p-plaquettes are shown in white (grey). For this tiling it 
is still true that e anyons can only be created and moved 
between neighboring s-plaquettes, and m anyons between 
neighboring p-plaquettes. Note that all strongly cou- 
pled plaquettes in this tiling are triangular. The Wp for 
these will therefore be three-body operators, making the 
code-ferromagnet coupling only a four-body term. On 
the other hand, weakly coupled plaquettes are octagons 
with eight-body Wp and nine-body terms required for 
the code-ferromagnet coupling. Since these many-body 
terms will most likely be generated by perturbative meth- 
ods, with a higher number of spins in a term generated 
by higher orders of perturbation theory, the difference in 
coupling strengths will arise naturally. 

Due to the practical difficultly in generating many- 
body terms, we can consider not coupling the octagonal 
plaquettes to the ferromagnet. Only the four-body terms 
required to coupled the triangles are then needed, which 
should be easier to implement than the five-body terms 
required for the square tiling. Despite the fact that only 
a fraction of the plaquettes are coupled to the ferromag- 
net, the memory is still stable against thermal errors. 
This is because any single spin error must still create at 
least one anyon on, or move anyons through, energeti- 
cally penalized triangular plaquettes. The energy barrier 
that increases linearly with system size is therefore still 
intact, and ensures that anyon creation and diffusion are 
exponentially suppressed. 

Unfortunately, stability against local Hamiltonian per- 
turbations does not remain strong without the coupling 
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of octagons. Without an energy penalty, two-body per- 
turbations are free to create and move anyons between 
next-to-neighboring octagonal plaquettes. This avoids 
the energy barrier and so leads to uncorrectable errors 
in a constant time. However, it is possible to avoid this 
by carefully considering what types of perturbation are 
present, and then designing the Wp such that they are 
unable to perform such hopping processes. For example, 
let us use Wp = -fp,i-fp^2-^p,3 triangular s-plaquettes. 
Here spin 1 is that shared with the neighboring tri- 
angular s-plaquette and the numbering proceeds clock- 
wise. Let us also use Wp = -fp,i-fp_2-^p,3 for triangular 
p-plaquettes with corresponding numbering. No near- 



est neighbor isotropic perturbation of the form I" I", for 
a S {x, y, z}, commutes with all of these operators. This 
means such perturbations will be suppressed by the en- 
ergy barrier and will not be able to move anyons between 
octagonal plaquettes. If only perturbations of this form 
are present in the system, the memory will remain stable. 
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